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Abstract: We probe the effectiveness and robustness of a simple gauge/gravity dual 

model of the QCD fireball that breaks conformal symmetry by constructing a family of 
I/") similar geometries that solve the scalar /gravity equations of motion. This family has two 

parameters, one of which is associated to the temperature. We calculate two quantities, 
^ the speed of sound and the jet-quenching parameter. We find the speed of sound to be 

universal and robust over all the geometries when appropriate units are used, while the 
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jet-quenching parameter varies significantly away from the conformal limit. We note that 
the overall structure of the jet-quenching depends strongly on whether the running scalar 
is the dilaton or not. We also discuss the variation of the scalar potential over our family 
of solutions, and truncate our results to where the associated error is small. 



1. Introduction 



The dual description of hot, strongly coupled gauge theories in terms of weakly coupled 
supergravity in a curved spacetime is one of the greatest surprises of the last dozen years. 
This correspondence is tantalizingly close to offering a theoretical framework by which 
one might hope to understand the dense matter produced at the Relativistic Heavy Ion 
Collider (RHIC) — after all, with collision energies around 40 TeV the produced plasma 
is certainly "hot" (on the order of a couple hundred MeV), and the quantifiable successes 
of hydrodynamic models applied to RHIC flow phenomena suggest that the matter may 
indeed be "strongly coupled". 

Unfortunately, it is not yet known what curved spacetime we should study in order 
to calculate in quantum chromo dynamics. The simplest examples of the gauge/gravity 
correspondence, like the one between M = 4 Super Yang- Mills and AdS§ x S" 5 , often 
exhibit decidedly un-QCD like characteristics. These include supersymmetry, exact con- 
formal symmetry and an absence of degrees of freedom transforming under the fundamental 
representation of the gauge group. 

In light of this somewhat bothersome obstacle, there remain several options for ap- 
plying these string theoretic techniques to hot QCD matter. First, one could continue 
calculating in AdS$ x S 5 and hope that the results obtained do not depend much on the 
details of the geometry. Such is the case for the ratio of the viscosity to entropy density 
obtained in [1] , in which the specifics of the geometry (and hence the corresponding gauge 
theory) "divide out" , leaving a result universal over a broad class of theories (for possible 
exceptions, see [2, 3, 4, 5]). Such results are probably rare, but highly desirable — their 
indifference to the details of the theory imply that if QCD has a gravitational dual, they 
hold there as well. However, since such universality can at best hold for some subset of all 
properties of hot QCD, it is useful to think about looking further. 

Alternatively, one may attempt to incorporate more QCD-like features into the corre- 
spondence (see for example [6, 7, 8] and references therein). The ideal situation, of course, 
is to find an exact solution to string theory or supergravity with more realistic features, 
such as non-conformality, less supersymmetry or fundamental matter. (The most ideal sit- 
uation would be to find an exact dual to QCD, but this seems at present too much too hope 
for.) A difficulty with this important path is that only a very small number of exact string 
or supergravity solutions at finite temperature with the appropriate features are known; 
for example, the M = 2* model [9, 10, 11] and the Sakai-Sugimoto model [12, 13, 14, 15]. 
Much important work has been done in this direction. 

Another possibility, however, is to sidestep the scarcity and complexity of exact solu- 
tions by simply postulating a spacetime geometry with the desired properties, preferably a 
simple one, without explicit knowledge of what matter content (if any) renders this back- 
ground a solution of the supergravity equations. For example, the metric of Kajantie, 
Tahkokallio and Yee [16] breaks conformal symmetry by introducing a simple dimensionful 
parameter in a warp factor while simultaneously turning on a scalar field. This metric 
was then explored by Liu, Rajagopal and Shi [17], and its consequences for a number of 
dynamical properties in the QCD fireball were worked out. If such models can yield useful 
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results, a great deal of effort obtaining exact solutions can be circumvented. 

Even with such a model, however, it is important to think carefully about the utility 
of results obtained via gauge/gravity duals. Although the model may move closer to 
having QCD-like properties, it is still not QCD, and whatever quantities one calculates 
are again most useful if they do not vary much over the class of spacetimes with whatever 
properties were imposed. For example, any result obtained for an arbitrary hot plasma 
lacking conformal symmetry that is not QCD is only useful for understanding QCD in as 
far as the result is not sensitive to which non-conformal plasma is considered. Put another 
way, "useful" results for non-conformal hot plasmas only need to be universal across the 
subset of gravity duals that break conformal symmetry. This idea is discussed at length in 
[17]. 

Thus one hopes to capture the properties of a whole class of gauge theories, ideally 
including QCD, by constructing a simple model with one or two desired properties. Two 
natural questions then arise. First, are any pathologies introduced into the gauge theory 
behavior by taking a gravity dual that is not known to solve any equations of motion? 
Even if one fully trusts the existence of the gauge/gravity correspondence, it only promises 
consistent and physical dual behavior for backgrounds that are genuinely consistent solu- 
tions to a theory of quantum gravity. It is possible that backgrounds that exist in the 
"gravity swampland" , without any extension to a consistent theory of gravity, would pre- 
dict unphysical behavior. And second, how robust are such models in calculating physical 
quantities of interest? 

To address these issues, in this paper we construct gravity solutions of varying temper- 
ature closely analogous to the solutions of [16] , with the same dimensionful parameter in the 
four-dimensional part of the metric, but which are known to solve the equations of motion 
for a single scalar field with a potential coupled to gravity [18, 19]; similar scalar /gravity 
models not designed to directly make contact with [16] include [18, 21, 22]. These solu- 
tions are necessarily more complicated than those of [16] . Moreover, these solutions come 
in a family parameterized not just by the temperature, but by an additional dimensionless 
quantity, here called a. We then study two properties of our metrics: a bulk property, 
the speed of sound, which we compare to the speed of sound calculated directly from the 
model [16], and the jet-quenching parameter q, which we compare to the results of [17]. 
Gratifyingly, we find broad agreement, suggesting that the use of the simpler toy model 
not known to satisfy equations of motion is not problematic in this case. 

Moreover, we are able to normalize the temperature scale such that the speed of sound 
for our family of solutions has the same dependence on temperature regardless of the 
additional dimensionless parameter a. However, we find that in this physical scale, the 
form of the jet-quenching parameter varies with temperature differently depending on the 
value of a. This suggests that our geometries may be quite robust in predicting a speed 
of sound; but on the other hand may not be as robust for calculating the jet-quenching 
parameter. 

The solutions we construct, although more complicated than those of [16], share the 
same warp factor for the four-dimensional part of the metric. It is not obvious in general 
that a whole class of solutions with different temperatures should exist with such a common 
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warp factor. Indeed, we find that the solutions we obtained for varying temperatures have 
slightly different values of the scalar potential. If the potential varied wildly between 
solutions, there would be no sense in which we could think of them as states of different 
temperature in the same theory. We show, however, that for the region of physical interest, 
the change in the potential is small compared to the potential itself, and thus the associated 
errors are small. This places a limit on the validity of our approach, where the variation 
of the potential becomes significant. 

Finally, in the model of [16], it was assumed that the one scalar field present was 
the dilaton, which affects the string metric perceived by the worldsheet in calculations of 
quantities such as the jet-quenching parameter q. A priori, the single-scalar solutions we 
construct are agnostic as to whether the scalar field is the dilaton or not. We compare 
the calculation of q between the two cases and find dramatically different behavior with 
temperature, in one case rising and in the other falling. We suggest that the evolution 
of the jet quenching parameter with temperature may provide a way of distinguishing 
whether the spacetime dual to QCD possesses a running dilaton, and hence potentially a 
strongly-coupled region somewhere in its extent. 

We construct our family of solutions in section 2. After reviewing the conformal 
symmetry-breaking model of [16] in section 2.1, we detail our solutions in section 2.2, and 
discuss the variation of the scalar potential in section 2.3. In section 3, we calculate the 
speed of sound both for the model of [16] and for our solutions, and show how with an 
appropriate choice of temperature scale, all models have the same functional form for the 
sound speed. In section 4 we calculate the jet-quenching parameter q for our models and 
compare to the result of [17], using both the temperature scale of [17] and the units where 
the speed of sound is universal. We also compare the behavior of q when the scalar field is 
treated as the dilaton to when it is not. We discuss our results and conclude in section 5. 

2. Probing the Plasma: Models 

In an effort to measure how useful a result computed in the gauge/gravity correspondence 
is to RHIC physics, the authors of [17] offer a pragmatic routine. First, choose a simple 
metric that preserves all the desired spacetime symmetries and reduces to A0IS5 in the 
UV, but breaks conformal invariance. Effectively, this allows one to "turn on" a feature 
of QCD not typically present in pure AdS/CFT. Next, one uses this metric as input to 
computation in the gauge/gravity correspondence. By studying how the results depend on 
the conformal symmetry breaking parameter, one hopes to measure the "robustness" of 
some result to the shift from model to QCD. 

In this section, we review the model of Kajantie, Tahkokallio and Yee (KTY) [16] 
employed by Liu, Rajagopal and Shi [17] to undertake this program, and then construct 
a class of backgrounds sharing the same conformal-symmetry breaking parameter that 
explicitly satisfy the equations of motion for one scalar coupled to gravity. 

2.1 KTY model 

Let us review the model of [16]. We use conventions where the action of the gravity /scalar 
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system reads 

5 =l6k/ d5 ^( fl -^* )2 - V "<*))' <21) 



In these conventions, the KTY model takes the form 



r: 



2 



Z 2 



ds 2 = i^/M 

with 



-h(z)dt 2 + dx 2 + 



dz 



2 1 



Hz) 



(2.2) 



4 



2f(z) = -cz 2 , h{z) = 1 - - A , (2.3) 

where the associated temperature is T = 1/ttzq. The model also contains a running scalar 
field $, 

$ = yi# 2 , (2.4) 

which was interpreted in [16] as the 5D dilaton; we will remain neutral on whether we treat 
this field as the dilaton or not until section 4. 1 If it is the dilaton, the string metric still 
has the form (2.2) but with 

2f str (z) = -cz 2 + = -cz 2 + . (2.5) 

The metric and scalar do not solve the equations of motion of the action (2.1), but it was 
supposed that some additional, otherwise-ignorable matter could be added to generate a 
true solution. 

There are three constant parameters associated to the solution: c, a parameter of 
nonconformality appearing in the warp factor; (j), the overall normalization of the running 
scalar; and zq, which controls the temperature. Only two are independent, however: we 
note that we can change coordinates: 

t — > Xt , x — > Ax , z — > Xz , (2-6) 

which is the isometry of pure AdS space (recovered in the limit c = 4> — > 0, zq — > oo) 
associated to scale transformations. The KTY metric is then unchanged if we also scale 

zq —> Xzq , c^X~ 2 c, 4>^X~ 2 4>, (2-7) 

corresponding to the freedom to rescale all mass parameters simultaneously. The two 
invariant dimensionless combinations can be taken to be, 

c4 « ^ . a = 4 • (2-8) 

In addition to the ratio of the temperature to the conformal symmetry-breaking parameter 
c, we also have the parameter a, encoding the ratio of conformal breaking in the metric to 



lr The canonical normalization for the 5D dilaton, used in [16], is $5 = ^/3/8$. 
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the breaking in the scalar Thus in addition to the temperature, there is a one-parameter 
family of conformal-symmetry-breaking geometries. The field theory interpretation of a is 
in general opaque. 

KTY fix a by appealing to the behavior in the meson spectrum of a different metric in- 
tended to correspond to the low-energy effective theory on the other side of the confinement 
phase transition [16]: 

20 , , 

&KTY = |g • (2.9) 

We shall consider this value, but we will generally leave a arbitrary, and examine the 
resulting physics for several different choices in the vicinity of (2.9). By considering this 
parameter, we will be able to quantify how much physical quantities vary as we go from 
one metric to another at fixed temperature, and hence get an idea of the robustness of our 
family of solutions. 

2.2 Model solving equations of motion 

We consider a metric ansatz with the general form, 

p2B(r) 

ds 2 = e 2A ^ (-h(r) dt 2 + dx 2 ) + ---dz 2 , (2.10) 

' n{r) 

where r is a suitable radial variable, along with a varying scalar $(r). In papers by Gubser 
et al., [18, 19] the equations of motion for this system following from the action (2.1) were 
worked out and found to be, 

A" -A'B' + - = 0, (2.11) 
6 

h" + (AA' — B')h! = , (2.12) 
2e 2B V + QA'ti + h(24A' 2 - 1) = , (2.13) 

in a gauge where the scalar field itself is used as the radial coordinate, r = $; we have 
neglected to write the scalar equation of motion which is not algebraically independent. 

In an attempt to generate a family of solutions to these equations with the same 
type of conformal symmetry breaking parameter as (2.2), we assume a warp factor for the 
four-dimensional directions of the same form as (2.3), which using (2.4) and the definition 
a = c/<j) becomes 

*ro = N(l/!4) < 2 -"> 

using $ itself as the radial coordinate. It is then straightforward [18] to solve (2.11) to find 
= log (?-\ + 1 ± 2 f log [ 1 + ) -log$--L*. (2.15) 



2 ) ' 2a 2 to ^ 'V 3 J to aVQ 

Note that integration constants in both A and B are undetermined by (2.11); we have 
fixed A to agree with the KTY metric, and B to reduce to AdS space with matching 



- 5 - 




Figure 1: Plots of the horizon function h(Q) for three different values of the temperature parameter 
hi, for a = 20/49. 



length scale R in the <1? —* limit. Beyond this parameter, this family of metrics depends 
on the dimensionless ratio a. 

We next see from equation (2.12) that the horizon function h(&) consistent with the 
equations of motion can be obtained by integrating [18] 

= ho- c 2 B?hi !^ d$'e- 4A ^ +B ^ , (2.16) 
Jo 

with the values of A($>) and B($>) obtained via (2.14) and (2.15). Here we chose the lower 
bound on the integral for convenience, since any change in it can be subsumed into the 
constants ho and hi, which are left undetermined by (2.12). We use the freedom in ho to 
impose h(<& = 0) = 1 to reproduce the AdS limit in the ultraviolet, so ho = 1. This leaves 
hi as a free parameter in the set of solutions, which we will associate with the temperature. 

For a black brane-type geometry we expect h to proceed monotonically from h = 1 at 
the boundary to h = 0, which defines the horizon; we are not interested in the space beyond 
the horizon. To verify that our proposed ansatz leads to actual black brane solutions, we 
must solve equation (2.16) for h(<&). As it turns out, for arbitrary a, the integral in (2.16) 
is quite complicated, but it can be evaluated for particular choices of a. The result for the 
KTY value of a = 20/49 can be obtained in terms of incomplete gamma functions, 



= 1 - Ahi 



„ /4801 801 \ „ / 4801 267(147 + 20^6$) \ . 

soor — , — -soor — , — ^ — '-} (2.17 

V 800 ' 800 7 I 800 ' 39200 J V ; 

„ /4001 801 \ „ /4001 267(147 + 20\/6$)\ 

801 r — , — + 801 r — , — i — '- 

V 800 800y ^ 800 39200 J 



where A is the enjoyable constant 



163840000000000 (§) 1/400 2 1 /i60 e 80i/800 



A = ^ ps 0.0017 . (2.18) 

264116234249604801 x 89V«» v ' 
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Figure 2: Temperature dependence on the horizon location. Deviations increase as T J yfc decreases. 



A simpler choice is the nearby a = 1/2, for which we find the polynomial result 

h($) = 1 - (l80$ 2 + 60\/6$ 3 + 45$ 4 + 2 V'u* 5 ) . (2.19) 

We studied a number of choices of a in the vicinity < a < 1, where the integrals could 
be solved; in the remaining sections we will also plot results for a = 1/5 and a = 4/5. 

For the particular value a = 20/49 we plot the horizon function for three different 
values of h\ in figure 1; we see that the function has the correct form, monotonically 
decreasing to the horizon, and the location of the horizon shrinks as h\ grows. Similar 
results hold for other values of a. Given that our solutions do indeed possess a horizon, 
we can use them to study a non-conformal plasma at finite temperature. 

We would like to translate the horizon location <§>ti into a value of the temperature. 
The temperature T of the geometries is related to the various metric functions by [18] 

T = — 1 v hn , (2.20) 

where the derivative acting on h is with respect to the horizon value note this is not 
the same as taking a derivative with respect to $ and then setting $ = since in general 
/i(<J>) is also a function of the parameter 

In figure 2 we show the relationship between the location of the horizon and the 
temperature for our solutions, and for the KTY model, for a = 20/49; very similar results 
hold for other a in the range 1/5 < a < 1. To evaluate this for the model of [17] is a simple 
exercise in differentiation. For the more complicated solution we study here, it requires a 
numerical routine, where the algorithm obtains the location of the horizon for many values 
of h\, finds the corresponding temperature using (2.20), and compares. The two agree in 
the large-T (small-c) region and deviate only slightly over the range of 1/2 < T/^/c < 1. 
This range was argued to be of greatest physical interest in [17], through comparison of 
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Figure 3: A plot of V(Q)/R 2 for three values of the temperature, and the corresponding horizons 
indicated as circles. The potential changes substantially, but most of the change is behind the 
horizon. 

the thermodynamics of the KTY model to lattice data. This small deviation is perhaps 
a first suggestion that switching to explicit solutions of the equations of motion does not 
create a substantial change. For smaller temperatures they diverge more substantially, for 
reasons we now describe. 

2.3 Variation of the scalar potential 

So far we have discussed equations (2.11) and (2.12), used to determine B(&) and /i(<I>) 
given A(&). However, we have not considered equation (2.13), which fixes V(3>) in terms 
of the other functions. As was noted in [18], in general this leads to a problem for solving 
the equations the way we have: there is no guarantee that any given set of functions A, 
B, and h varying with temperature will lead to the same scalar potential V(<&). If the 
potential changes, then the solutions cannot be thought of as belonging to the same gauge 
theory dual; instead, as one varies the temperature, one is also changing parameters of the 
Lagrangian, and producing a vacuum of a different temperature in a different theory. 

The solutions to the equations of motion we find, assuming the KTY ansatz for A{§) 
and using (2.13), indeed do not lead to the same potential for different temperatures. While 
it is possible that there exists some arrangement of matter fields for which the potential 
remains fixed, it is in principle difficult to find. We could also consider a more complicated 
ansatz, where we do not assume a simple form for -A(3>), but in addition to the issue of 
complexity, this defeats the purpose of comparing to the KTY model. 2 

Here we take a practical approach, and see how much it matters. Examining how much 
the potential is varying with the temperature (that is, the location of the horizon), we find 
the following. As T changes, V(0) = —12B? remains fixed, while the largest variation 
occurs at larger However, points at sufficiently large are always behind the horizon 

2 Solutions to the gravity /scalar system with fixed scalar potential have been studied recently in [21, 22]. 
These solutions indeed have different warp factors for different values of the temperature. 
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at 5>h, and hence the region of largest variation does not contribute to the physics. As 
the temperature increases, the region where V(<&) is varying substantially retreats towards 
$ = 0, but at the same time <&h retreats in front of it, tending to keep the region of strong 
variation harmlessly behind the horizon. 

We can quantify the change in V{<&) at some point <I> as we vary the solutions over 
the widest physically relevant range: from T = to the temperature where that point 
disappears behind the horizon. We find that the variation in V at a given value of the 
scalar increases monotonically with $; thus for the solution at a given temperature, the 
point that varies the most over all temperatures to which it contributes is always the point 
with the largest value of 3>, namely &h itself. We can thus parameterize the variation of 
V for a solution at a particular temperature T, by how much the potential at the horizon 
point V(&h) varies between that solution and the zero-temperature solution. 

We have performed this analysis, and for a = 20/49 for example we find the values 



AV(® h ) 



0.005, 



T/y/c=l 



0.04, 



(2.21) 



T/y/c=l/2 



over the range in T '/ ' y/c argued in [17] to be most physically relevant. Thus although 
the solutions to the equations of motion we find are not solutions for precisely the same 
potential, they are sufficiently similar to lead to only a small error of a couple percent over 
the region we are interested in. This error, moreover, is over the entire possible range of 
temperatures; for a practical calculation comparing two solutions with similar temperatures 
the error will be smaller. 

In general we find that smaller values of T j y/c have worse behavior; thus the larger the 
conformal symmetry-breaking, the more the effect matters. Similar results hold for other 
values of a, and as a increases, the "safe" region moves further in towards smaller values 
of T j y/c (and vice versa). In the sections that follow, we arbitrarily choose a largest error 
in AV/V that we will tolerate of five percent, and truncate the results for various values 
of a each at this cut-off. The remaining points should be robust to the variation of the 
potential to within a few percent. 

This is a trade-off for attempting to find exact solutions that, while more complicated 
than the KTY model, are still relatively simple; in principle there should exist solutions 
with precisely the same value of V(<3?) for all temperatures, but these are in general hard 
to obtain, and additionally, most of these solutions will bear no particular resemblance to 
the KTY metric for any value of c. 

It is also possible to vary our one other parameter, a, as we vary T/y/c, to try and 
mitigate the variation of the potential; that is, we can try to find the path through the 
two-dimensional parameter space of T/y/c and a that minimizes the variation of V($>). 
It is possible to decrease the variation of the potential by several orders of magnitude by 
making a compensating change in a, though one cannot cancel it precisely. We will not 
attempt to systematically perform this compensation, since the error from simply changing 
T/y/c is already small for the region we are most interested in. 
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3. The speed of sound 



Having our solutions in hand, we would like to compare their predictions to those of the 
KTY model, to try and see whether the use of metrics not solving known equations of 
motion still leads to reasonable and robust results. We will compare two quantities: first, 
a "bulk" quantity not involving the string worldsheet, the speed of sound. In the next 
section, we shall consider a worldsheet-related number, the jet-quenching parameter. 

In a non-conformal plasma, one expects richer thermodynamic behavior than in the 
conformal case. For example, conformal invariance fixes the speed of sound, defined in 
terms of the pressure P and energy density e as 

* = (3-D 

at (? s = 1/3. This condition comes from the tracelessness of the energy-momentum tensor, 
as mandated by invariance under dilations. In a non-conformal plasma, on the other hand, 
c 2 s can vary. 

A convenient formula for the speed of sound is 

2 _ dlogT _ dlogT f dlogs y 1 
° s dlogs' d<S> h \ d* h J ' [6 - Z) 

valid in the limit of zero net baryon density, a reasonable approximation in the heavy-ion 
fireball. Besides the temperature it requires knowledge of the entropy density s, which for 
our systems is simply proportional to the horizon area of the black brane, taking the form 

e 3A($ h ) 

S =W < 33 > 

and since both T and s vary with it is easiest for us to compute (3.2) in terms of the 
variation of each quantity with <£/j, as indicated. 

For the KTY model, the speed of sound is straightforward to evaluate: 

C -3TTW < 34 > 

Note that this formula is independent of a, which for the KTY model entered only into 
the scalar field and not into the metric (which controls s and T). The result is plotted in 
figure 4. 3 This expression clearly has the right limit of c 2 s — ► 1/3 as we restore conformal 
symmetry c — > 0, and it monotonically decreases as the temperature gets small to approach 
zero in the limit of zero T. This decrease below the conformal value is reasonable on physical 
grounds. 4 



3 A very similar functional form for the sound speed was found in [23], based on the model in [24] which 
resembles [16] without a running scalar. 

4 For example, in the simpler case of the perfect fluid P — toe, we have c 2 a = w and the interpolation 
between c? s = 1/3 and c? s = is just the progression from ultrarelativistic to non-relativistic species. This 
plasma would have no quasiparticle description, but the behavior of the speed of sound is analogous. 
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Figure 4: Behavior of the squared speed of sound with respect to the dimcnsionlcss quantity T/y/c 
plotted for several values of the parameter a. Also shown is the conformal value (c 2 s — 1/3), and 
the result given by the KTY model. 



For the family of solutions obtained in this paper, calculating the speed of sound in 
general requires a numerical routine. The algorithm is simply to scan over many values of 
the horizon location by varying hi, and calculating the resulting values of the temper- 
ature and entropy. We then compare neighboring values of $/i to calculate the speed of 
sound. Note that because this calculation involves comparing two solutions with distinct 
temperature, it is particularly sensitive to the variation of the potential between those so- 
lutions; we should not expect this to be a reliable calculation in the limit that the potential 
variation becomes substantial. 

Crucially, the results of this calculation vary depending on the value of a chosen. In 
figure 4 we plot the results for a few representative values of this parameter, along with the 
single result (3.4) for the KTY model. For our results, as each model approaches smaller 
values of T/y/c the error due to the variation of the potential becomes larger. We have 
therefore imposed a somewhat arbitrary cutoff, including only results with AV/V < 0.05. 
Following the results beyond this cutoff eventually results in pathological behavior, in 
particular the speed of sound squared going below zero at nonzero temperature. 

For all values of a, the speed of sound asymptotes to the conformal value at large 
temperature and decreases at smaller temperature, analogous to the KTY model. The 
curvatures of our solutions are slightly different from KTY, being initially natter and 
then dropping more abruptly as the temperature is lowered. Smaller values of a give 
a lower speed of sound at fixed temperature, while larger values give a greater speed. 
Unfortunately, the variation of the potential and the associated curve truncation do not 
allow us to determine the small T/y/c behavior of the speed of sound. As discussed in 
[25], lattice simulations for hot QCD matter predict a "strong" drop in the speed of sound 
near the critical temperature (T c ~ 170 MeV), followed by a rise below it. Obviously, 
such behavior is not present in the KTY model, which monotonically approaches c 2 s = 
as T/y/c — > 0. This is because we work exclusively with the "deconfined phase" geometry, 
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Figure 5: Speed of sound curves measured in rescaled units T/y/c^. In these units, our model 
makes one prediction for the dependence of the speed of sound on the ratio of temperature to a 
dimcnsionful combination of the parameters <j> and c. We also compare to the KTY result and an 
analogous function depending on T 4 . 



which corresponds to T > T c . For an attempt to realize this behavior in a gravity dual, 
see [18, 20]. 

We have plotted our results in terms of T/y/c in order to make comparisons with [17]; 
they naturally used c as a reference parameter since only it appears in the metric, which 
controls the quantities of interest. For full solutions to the equations of motion, however, 
the parameter <fi feeds back into the metric as well via the parameter a. It is not obvious 
that T/y/c, as opposed to T/y/fy or any combination of the two massive parameters, is the 
most useful scale to employ. In other words, our additional dimensionless parameter makes 
it unclear which units are the natural ones to use. 

We can use our results for the speed of sound to try and resolve this issue physically. 
The speed of sound is a calculable quantity with physical consequences. We may choose, if 
we can, to set our scale for each model precisely so that in these units, the same temperature 
always gives the same value of the speed of sound. That is, we can seek to measure T 
relative to whatever combination of c and <j> we need to make the speed of sound curve look 
identical for all values of a. 

We find that plotting the speed of sound as a function of T/y/c^, where c a is an in- 
dependent rescaling of c given by (c 4 / 5 , Ci/ 2 , c^o/aq, c i/5) = (0.14c, 0.32c, 0.46c, 1.05c) with 
similar scalings for other values of a, creates a convergence between the various curves of 
our model, all agreeing substantially with each other over their region of validity. This 
suggests that the appropriate units to measure our results with respect to are T/y/c^, 




which we will do. Since this rescaling depends on a, we can think of c a as a combination 
of the massive parameters c and <fi. In figure 5, we plot our rescaled sound speeds, along 
with the KTY result (3.4) and also the function 




which is a better fit to our unified result for the speed of sound. 

We can compare our result to the calculation of the speed of sound in other models, 
such as that for the M = 2* theory in [26]. In terms of conformal symmetry breaking 
parameters nif and m^, they found to leading order in each mj/T and m^/T, 

-aO-W?)"-!^?)'*-")- 

which broadly matches KTY and our results in the leading behavior: all decrease from 
the conformal value as the temperature goes down. For generic values of mj and m& 
(including the M = 2 limit mf = nib), the (nib/T) term will be subleading, and in this 
case the M = 2* speed of sound matches the KTY model to order 1/T 2 (including the 
sign) as long as we identify the symmetry-breaking parameters c = 2[r(3/4)] 2 m 2 /37r 2 . In 
the limit m/ <C nib, however (including m/ = 0), the N = 2* models predict the first 
subleading term to be proportional to —1/T 4 , which is a better match to our models and 
the form (3.5). It is possible that the best fit to our models contains both terms at leading 
order, as in (3. 6). 5 

We have seen that as long as we use the right units, our determination of the speed of 
sound is quite robust, and broadly similar to the KTY result though not of precisely the 
same form. As a consequence, we can use the physical quantity of the speed of sound to 
circumvent our ignorance about the physical meanings of c and 4> and use the running of 
the speed of sound itself as our effective scale. Let us now examine a worldsheet-dependent 
quantity, the jet-quenching parameter q. 

4. Jet-quenching parameter 

The authors of [17] calculate a number of observables associated to the propagation of 
quarks in the plasma defined by the KTY model, which on the gravity side are associated 
to worldsheet calculations. Here we subject our class of models to the calculation of one 
of these, the jet-quenching parameter q. 

The parameter is computed via the integral 6 

1 / Z"** d$ 

q ~ na' yj y/h($)[l - exp[6A($) - 25($)] 

where a' is the string tension, and the bars above the metric functions are reminders that 

they must be transformed into string frame. We can calculate this quantity for our model 

as well. Note that to do this, we must make a choice as to whether $ corresponds to the 

dilaton or not, since (4.1) is sensitive to the string metric. Seeking to have things all ways, 

we will try both choices and compare. 

5 Other theories that approach (? a = 1/3 from below in the high-T limit include the model of [23, 24] which 
is closely related to KTY, the finite-temperature Klebanov-Strassler cascade [27, 28], the Gubser-Nellore 
model [18], and the Giirsoy-Kiritsis-Mazzanti-Nitti solutions [21, 22]. Models with constant c? s < 1/3 also 
exist, including the Sakai-Sugimoto model [29], the Chamblin-Reall models [18] and the model of Springer 
[30, 31]. 

6 The definition used by [17, 34, 35], as proposed by [32, 33], has been called into question by [36]; since 
our primary motivation is to compare to the [17] computation we will use their definition. 
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Figure 6: Jet-quenching parameter plotted against c/T 2 for several values of a, and in the KTY 
model. The definition of the quantity A follows [17]. 



First we calculate q for various values of a assuming the scalar is the dilaton, and 
plot it with respect to T '/ ' ^fc in figure 6, along with the value obtained by [17]. Following 
[17], we show q relative to T 3 V% taking into account that many energy loss models used 
to describe jet quenching at RHIC take q to scale like T 3 [37]. We notice a few things: 
first of all, the high-temperature (conformal) limit is robust, with the KTY model and our 
solutions for various a all approaching the same value at q/(T 3 \/~\) ~ 7.5. However, the 
functional forms of the various values of a diverge substantially as conformal symmetry 
breaking is introduced, with small a possessing the greatest slope. 

We established in the previous subsection, however, that a more physical scale is to 
measure quantities relative to T/^/c^, since in this scale, every value of a has the same 
speed of sound. Thus if the jet-quenching parameter varies when measured relative to 
T/^Ca, this represents a physical difference between models and not just an unknown and 
arbitrary rescaling of the reference scale. 

In figure 7 we plot the same results relative to the physical scale. Interestingly, we still 
find distinct curves for q at different values of a. Thus the variation of a represents a true 
physical deformation, allowing us to vary the functional form of the jet quenching while 
keeping the "bulk" property of the speed of sound fixed. Notably, the relative positions of 
the different values of a have flipped, with large a now having greater slope. 

We can compare this result to the case where the scalar is not treated as the dilaton, 
so the string metric warp factors are identical to the Einstein cases (2.14), (2.15); the results 
are plotted in figures 8, 9. We note that the functional form of q changes substantially 
between the two cases. For a running dilaton, regardless of a the ratio g/(T 3 \/A) increased 
as the temperature dropped; for the absence of a running dilaton we find the opposite, 
with the ratio generally falling, more pronouncedly in the scaled (physical) coordinates. 
Again the switch to scaled coordinates T '/ \fc^ also exchanges the ordering of the models. 
We discuss these results in the final section. 
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Figure 7: Jet-quenching parameter in scaled coordinates. The universality observed in the speed 
of sound calculation is no longer present; our model does not seem to uniquely determine the 
temperature dependence of q. 
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Figure 8: Jet-quenching parameter when 4> is not the dilaton, in unsealed coordinates. 



5. Discussion and conclusions 



We set out to try and understand whether the use of simple model geometries not known 
to solve any equation of motion but possessing a key desired characteristic was reasonable 
and robust. At this point we can attempt to address these questions. 

There is no indication whatsoever that the use of the KTY model leads to unphysical 
results, or indeed results particularly at odds with the class of solutions to the scalar /gravity 
system we construct. For both the speed of sound and the jet-quenching parameter, the 
KTY model makes predictions solidly in the middle of the range of results of our solutions, 
and with similar functional forms. So this is reassuring: if there is a hidden difficulty with 
such a simple model, we have not found it. 
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Figure 9: Jet-quenching parameter when $ is not the dilaton, in scaled coordinates. 



Moreover, we have found that the functional form of the speed of sound is universal 
across our models, regardless of our tunable parameter a, as long as we measure it in the 
right units; it also agrees with the leading-order behavior of a class of M = 2* models [26] . 
Thus the predictions of the sound speed are independent of whichever representative of 
our solutions one chooses to use, and hence can be said to be significantly robust. This 
does not imply that they are universal across all non-conformal plasmas, but are part of 
a broader class including the KTY model and others that approaches the conformal value 
from below as the temperature gets large. 

For examining the results for the jet-quenching parameter q, it is useful for us to 
distinguish two different notions of "robustness" . As a first point of comparison, we can 
look at the q results in the context of the results of Liu, Rajagopal and Yee [17]. In this 
work, a dominant thrust was to identify whether or not a given result was "robust" to the 
introduction of non-conformality. The idea is that a "robust" quantity would change its 
value little as the conformal symmetry breaking parameter relative to the temperature was 
varied across some physically interesting range. Although the variation of the potential 
does not allow us to scan an arbitrary range of, say, c/T 2 , we can nevertheless make an 
comparison of most of our models over the range considered in [17] . 

For the a values 20/49, 1/2 and 4/5, we can follow q all the way from c/T 2 = to 
4. Over this range, one finds that the value of q/T 3 y/\ changes by about 15% for all a 
values. By comparison, over the same range the KTY result was found in [17] to increase 
by about 28%. From figure 6 it is easy to see that the largest variation is expected from 
a = 1/5. Unfortunately the variation of the potential limits us in this case to the range 



< c/T 2 < 1.4, across which q changes by 



o. 



When the scalar field is not the dilaton, 



figure 8, the jet quenching parameter will again change most when a = 1/5 — here 13 % 
for < c/T 2 < 1.4. When a = 4/5,1/2 or 20/49, we can explore the range < c/T 2 < 4, 
for which we find a 5%, 10%, and 10% change respectively. Accordingly, over these ranges 
in c/T 2 , the value of q for a given a is relatively robust to the introduction of the conformal 
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symmetry breaking parameter c following the definition of [17]. 

There is another important way in which our results may be judged for "robustness" . 
It is evident that none of the geometries we have studied is QCD itself, so their predictions 
are valuable exactly as far as they are not particularly dependent on the details of the 
model picked, as we have discussed. The ideal would be to pick a simple model with a 
few constraints — no conformal symmetry, the right kind of matter, and so on — and 
hope that the physical results (for at least some quantities) would then be independent of 
details of which model we have picked — in our case, of which value of the dimensionless 
parameter a was used. We found that the choice of a could be used to fix a unit system 
natural for the bulk property of the speed of sound, so we can normalize the variation of the 
sound curve across all models. Even with this normalization, however, the jet-quenching 
parameter had substantial variation even over the narrow window of a-values we were able 
to explore. 

Thus, on the one hand we find that the speed of sound is quite resilient. All of our 
models lead to the same functional form for the sound speed over their regions of validity, 
which can be made to coincide in the proper set of units; they closely resemble the model 
of KTY [16] as well. 

On the other hand, although our calculation of the jet-quenching parameter seems quite 
physical, it is less robust in this second sense, at least away from the conformal limit, as its 
functional form varies significantly with a. Without knowing the true field-theory dual, it 
is impossible to know what a corresponds to, and hence which value should be picked, and 
therefore which curve to use. Moreover, we could have chosen a more complicated starting 
point for our ansatz with more dialable parameters, and there is no reason to think the 
variation of q would not also persist over such a multi-dimensional parameter space. 

Thus if we want to extract universal properties of something like the jet-quenching 
parameter, we either have to provide new input constraints into our models, or look for a 
property insensitive to fine details, caring only about a gross behavior of the quantity. One 
example of the latter might be the substantial difference in the behavior of q between when 
<]? was treated as the dilaton and not. Although the details of the curves vary, one could 
try to extract an overall lesson that a running dilaton leads to a growth of g/T 3 \/A as the 
temperature diminishes, while a constant dilaton leads to the opposite. Currently, many 
energy loss models used to describe jet quenching at RHIC take q to scale like T 3 [37], 
giving no preference for the fate of our scalar field. Determining higher order temperature 
dependence would constitute another input into future models of the plasma, and place 
constraints on the available freedoms. 

In conclusion, we have found no reason to avoid simple models of the QCD plasma, 
but we have reinforced that the robustness of their predictions must be handled with care. 
Some quantities such as the sound speed may be quite robust, while others such as the 
jet-quenching parameter may be somewhat robust over the variation of a single parameter 
but less so once additional dialable parameters are introduced to the family of models. 
Even so, they may contain useful information such as the relationship between the dilaton 
and the global behavior of the function q(T). Further study of these questions in the future 
is surely warranted. 
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